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ABSTRACT 


In every wake of life, the flow of fluids through pipes is encountered. The major problem encountered while analysing pipe flow problems is obtaining friction factor. 
Though Moody's diagram helps evaluate the friction factor, the obtained solution is error-prone due to errors in reading the graph. So, to remove the mistakes using 
hand calculations and improper use of diagram, an attempt has been made in this research article to automate the process of pipe flow modelling. The Colebrook-White 
equation has been iteratively solved to obtain the friction factor. The modelling is done using Python as it is easy to use and has a vast library backup. The robustness of 
the developed program has been demonstrated by plotting Moody's diagram using the code. Three different pipe flow problems through a single pipeline are solved 
using the developed program, and the obtained results are precisely equal to those shown in the literature. 
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Nomenclature: 
A Cross-sectional area of pipe (m’) 
c, Fanning's friction factor 
d_ Pipe diameter (m) 
D,, Hydraulic diameter (m) 
L_ Pipe length (m) 
QO Volume flow rate (m*/s) 
R, Reynolds number 
v_ Average flow velocity 
f Darcy's Friction factor 
p Fluid density (kg/m’) 
Ap Pressure drop (Pa) 
e€ Wallroughness(m) 


1. INTRODUCTION: 

Understanding the fluid flow regime is of utmost importance from the engineer- 
ing point of view when a fluid flows. Reynold's number (Re) [1], which is the 
ratio of inertia to viscous force, talks about the flow regime, i.e., whether it is lami- 
nar or turbulent [2]. For Pipe flow, if the Re is less than 2000, it is called laminar, 
whereas when it is more than 4000, it is called turbulent flow [3]. Laminar flow is 
where laminas of fluid glide smoothly one over the other. In contrast, the turbu- 
lent flow is the one intermingling of fluid laminas. In other words, turbulent flow 
can be said to be where the flow becomes an irregular function of time [4]. And as 
the flow velocities in turbulent flows become periodic with respect to the time, 
we cannot have a steady flow. On the other hand, these fluctuations are random, 
so a statistical description is possible. The variation may be brought down by the 
law of probability so that an average velocity is defined [5]. 


Usually, all flows in nature are turbulent, so laminar flow occurs with minimal 
velocities, which are not usually encountered in practical applications. So the 
most vital thin in turbulent flow is that the fluid friction viz. layer to layer and 
between fluid and sold friction gets enhanced [6]. Which means it offers more 
resistance for the given type of fluid. This means for a given flow rate through a 
duct, the pressure drop across a length of a duct is much more than that of a lami- 
nar flow. Fora fluid flow in a pipe (Figure 1), the wall shear stress (t,,) can be eval- 
uated as [7]: 


Figure 1: Force balance of a fluid element in a pipe flow 


This means that stress is responsible for the pressure drop in the pipe. Now, if the 
above expression is replaced in the skin friction coefficient (c,), then it is also 
called the Fannings friction factor [7], as shown in Eq (2). 


Where v is the average flow velocity (v=Q/A), Darcy defined his friction factor 
as four times that of Fanning’s viz [8]. 


Therefore, the pressure drops (Ap) across a pipe length L in a fluid of density p 
and velocity v is given by: 


And the power to be supplied to the fluid (of flow rate Q) to overcome is pressure 
drop is given by: 


(5) 


Power = Ap XQ 


If the flow is laminar, one can quickly get the relation between friction factor f 
and Re, as shown in Eq. (6). 


But when the flow becomes turbulent, friction between the fluid and solid wall 
and within the fluid enhances due to mixing in the turbulent zone. So due to arise 
in flow resistance, the pressure drop gets enhanced. But there is not a single equa- 
tion by which one can fit the pressure drop and flow rate relation or f-Re relations 
through the entire range of turbulent flow regimes. So, to handle such situations, 
experiments come in. People like Stanton [9], Nikuradse [9,10], etc., have 
worked on developing such expressions, and finally, Moody documented them 
in the form of a diagram which is famously known as Moody's diagram [11]. 
They have observed that in turbulent flow, the fis not the only function of Re but 
wall roughness as well. But reading this diagram requires a great deal of interpo- 
lation, which can result in an error. So, to help the researchers working in the field 
of pipe flow, some sort of computations is required, which will save time and 
make the error of the computation free. 


Python is an easy programming language with a light and user-friendly syntax 
[12-15]. Moreover, when it comes to numerical computations, its modules viz. 
NumPy and SymPy are of great help [16-20]. NumPy has a multi-dimensional 
array and matrix data structure and has compact storage. Moreover, when it 
comes to loops, the NumPy is swift. Moreover, the slicing aspect of the array is 
compelling when it comes to numerical computations. Matplotlib.pylab [21] isa 
convenient library for data plotting. Pylab also uses NumPy, so while using 
Pylab, calling NumPy is not required. 
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In this research article, a python-based approach has been used to solve the fric- 
tion factor and the development of Moody's diagram. Also, three fundamental 
pipe flow problems are solved using the functions developed in Python. 


2. COLEBROOK-WHITE EQUATION AND THE EVALUATION OF 
FRICTION FACTOR: 

For fully developed turbulent flow in a pipe, the Darcy friction factor (f) can be 
approximated with the help of the Colebrook-White Equation [22], which is as 
follows: 


Where D, is the hydraulic diameter which gets reduced to the pipe diameter d for 
the circular pipe. Eq. (7) is an implicit equation which cannot be directly solved 
for the f So, an iterative procedure can be adapted to solve for/ The steps to solve 
the Eq. (7) implicitly are as follows: 


_ 


Rearrange the Colebrook equation into the following form: 


2 
= Dh ae (8) 
€ 2.51 
-2log 354 mH) 


ti. Choose the initial guess for /(say /,). For a better guess, use a value between 
0.1 and 0.01, as this is the range of f. 


iii. Solve Eq. (8) fornew/(say/,) by using /, as: 
2 


1 


fr = 
-21og € 2.51 
(55 7d* elfg 


iv. Obtain the error (|/,-/, |)and check for convergence. 


v. Update the /, by obtained value /, as anew guess for the next iteration. 
vi Repeat steps iii to v till the convergence is achieved. 


3. IMPLEMENTATION OF COLEBROOK EQUATION VIA PYTHON: 
The algorithm discussed in the above section has been implemented in Python 
via following function: 


def £ f£(e,d,Re,f_g): 
if Re<2000: 
£ n=64/Re 
else: 
error=1 
count=0 
e_d=e/d 
while error>1.E-5: 
count+=1 


a=e c_a/ (3 ) +2. 51/ (R 
= n=(1/ (- pans: nen 


f Error caiculati 


error=abs (f n-fg 


return fn 


First, Re is checked; if it is less than 2000 then Eq. (6) is used to evaluate the /, 
whereas the Colebrook equation is solved for other values. Moreover, one can 
also observe that in the equation, relative roughness is used, which is the ratio of 
wall roughness to that of pipe diameter (€/d ). Therefore, to use this function, one 
has to supply the wall roughness, pipe diameter, Re, and guess value for friction 
factor as shown in the below-mentioned code snippet. 
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d=150*10** (- 3) 


Re= 5000 


print (f! f= [round (£ric fact, ol) 


The output of the above code is f= 0.038496, which is the exact value of the fric- 
tion factor for the given arguments. 


One can also go one step further and plot the value of f for different Re and é/d. 
Below code will produce the result, which is nothing, but Moody's diagram and 
Figure 2 shows the diagram obtained from the code. One can observe from the 
plot that for every roughness in the turbulent region, there exists a Re beyond 
which the dependence of f on Re ends, i.e., the line becomes horizontal. 


from pylab import ~* 


# Range of a : ness 


e= Sian eesiiy E-8,0. 004, 10) 
da=0.05 


of Reyno 3 k 


R= =arange (700, 1. _E7, 100) 


figure (2,dpi=300) 
for eps in e: 
fric_fact=empty(len(R)) 
£ g=1 
for i,Re in enumerate (R): 
if Re==2000: 
fric_fact[i]=100 # Some arbitrary value (ver 
else: 
fric_fact[i]=f_f(eps,d,Re,f_g) 


semilogx(R,fric_ fact, '—', label=f£"e/d-{eps}') 
if eps/d<0.1: 
text (R[-1],fric_fact[-1],f" 
text (R[-1],¢[-1]/d+0.015,' e/d') 


{round (eps/d,7)}') 


xlabel ('Re', fontsize=14) 
ylabel ("£', fontsize=14) 
aieate poe; 0.1) 
$gqrid( 
for pos in tf right", ‘top"]: 
gca() . spines [pos] .set_visible (False) 


savefig ("Moody's.jpg") 
show () 
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Figure 2: Friction factor vs Re for different values of relative 
roughness. 


4. CLASSES OF PROBLEMS TYPICALLY ENCOUNTERED IN A 
FLUID FLOW INASINGLE PIPELINE: 

When flow through a single pipeline happens then, three distinct types of prob- 
lems usually are arisen, which are listed below: 


I. The pipe diameter, pipe length, and flow rates are given, and one has to find 
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the pressure drop and corresponding pumping power required. 


ii. Head loss over a pipe length for a pipe of known diameter is given; the target 
is to find the flow rate and power required. 


iii. The flow rate and the head loss over a given pipe length are known the pipe 
diameter has to be evaluated. 


The first class of problem is explicit, and one can quickly solve the pressure drop 
and power using Eq. (4) and (5). However, the second and third types of prob- 
lems cannot be solved explicitly as velocity and friction factors are unknown in 
the second, and diameter, and friction factors are unknown in the third class. So, 
these problems can only be solved by iteration procedure. To solve these prob- 
lems with the help of Python, following apart for the function to evaluate friction 
factor, the following more functions are created to ease up the problem solving: 


This function evaluates Ah, Ap 
and accept f,L,d p,v,g as 
input arguments. 


def press drop(f,L,d,9,v,g) : 
dp=£* (L/d) * (p*v**2) /2 
dh=dp/ (p*g) 


return dp,dh 

This function evaluates power | def Power (a,b): 

with Ap and Q as input return a*b 

This function retums flow | def vel _Q(Q,d): 

velocity when discharge and return Q/ ( (pi/4) *d**2) 
diameter are given 
This function retums Re for | def Rel _num(p,v,d,p): 

given p,v,d, u return p*v*d/p 

This function returns Re for| def Rel num v(v,d,v) =: 

given v,d, v return v*d/v 

This function returns velocity | def v_Ah(Ah,d,L,g,f_g): 

when Ah is given f are known return sqrt (Ah* (d/L) *2*g/f_g) 

This function retums diameter | def d_Q Ah(Q,Ah,L,g,f_g): 

when Q and Ah are known return 

((£_g/Ah) * (L/ (2*g) ) * (4/pi) **2*Q**2) ** (1/5) 


Now let us demonstrate how these problems can be solved using the above 
python functions. 


i. Determine the frictional head loss and pumping power required when water 
flows through a 300 m long steel pipe of 150 mm diameter at a flow rate of 
0.05 m’/s. The kinematic viscosity of water and the wall roughness of the 
pipe are v= 1.14x10° mm ande=0.15 mm [23]. 


Python Solution: 


d=150*10** (-3) 
O=0.05 
v=1.14*10** (-6) 
€=0.15*10** (-3) 
o=1.0*10**3 
p=p*v 


v=vel_Q(Q,d) 


Re=Rel_num(p,v,d, ) 


friction_factor=f_f(e,d,Re,0.1) 


tnrinti'ft = '.Ff 


print rt 


Ap, Ah=press_drop(friction_factor,L,d,p,v,g) 
print ('Ah = ',round(Ah,3), 'm') 


p=Power (Ap, Q) 
print ('Power = ',round(p,3), 'W') 


The output of the program is: 4h= 16.67 mand Power= 8176.72 W 


ii. Oil ofkinematic viscosity of 10° m’/s flows ata steady state through a pipe of 
diameter 100 mm and of surface roughness of 0.25 mm. Fora pipe length of 
120 m, the friction head loss is 5 m. What is the flow rate of oil through the 


pipe [23]. 


International Education & Research Journal [IERJ] 


Python solution: 


v=10** (-5) 
d=100*10** (-3) 
€=0 .25*10** (-3) 
L=120.0 

Ah=5 


count=1 

while error > 1.E-5: 
v=sqrt (Ah* (d/L) *2*g/f_g) 
Re=Rel_ num _v(v,d,v) 
f newFf f(e,d,Re,f_g) 


error=abs (f_new-f_g) 
£ g=f_new 


count=count+1 


v=sqrt (Ah* (d/L) *2*g/f_new) 
O= (pi/A) *d**2*y 

print ('iteration = 
print ("Discharge 


', count) 
",Q) 


The program output is: Q= 0.0126 m’/s. 


iii. Determine the size of a pipe (c=0.15 mm) needed to transmit water 
(v=1.14«10° m’/s) to a distance of 180 m at 0.85 m’/s witha loss of head of 9 
m [23]. 


Python solution: 


{ 


L=180.0 
O=0.085 
Ah=9.0 
v=1.14*10** (-6) 
€=0.15*10** (-3) 
g=9.81 


£ g-0.1 


error=1 

count=0 

while error>1.E-5: 
d=d_Q Ah(Q,Ah,L,g,f_g) 
v=vel_0(Q,d) 
Re=Rel_num_v(v,d,v) 
£ new=f_f(e,d,Re,f_g) 
error=abs (f_g-f_new) 


£_g=f_new 
count=count+1 
print (count, error) 


print ('final diameter =',d) 


The output of the program is: d=0.1872987 m. 


It has been observed that the results obtained by the computer programs are the 
same as those obtained in the literature [23]. 
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5. CONCLUSION: 

In this manuscript, pipe flow modelling has been done using Python. A detailed 
explanation of the pipe friction factor has been done. A precise algorithm has 
been explained to obtain the friction factor in the turbulence zone using the Cole- 
brook- White Equation. Moreover, three classes of fluid flow problems ina single 
pipeline were presented and solved using the functions developed in Python. The 
capacity of function to obtain friction factor has also been used to plot Moody’s 
diagram. The methodology provided in this article will help practice engineers 
and researchers solve pipe flow problems without errors. The developed com- 
puter programs are very robust and can be extended to solve any flow problem. 
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